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Abstract 

D , We compute the renormalized four-point coupling in the 2d Ising 

' model using transfer- matrix techniques. We greatly reduce the sys- 

tematic uncertainties which usually affect this type of calculations by 
^ \ using the exact knowledge of several terms in the scaling function of 

?-h ' the free energy. Our final result is g± = 14.69735(3). 



1 Introduction 



An important information on the physical properties of a quantum field the- 
ory is given by the renormalized four-point coupling, which is defined in terms 
of the zero momentum projection of the truncated 4-point correlator. At the 
same time, if one is interested in the lattice discretization of the theory, 
this renormalized coupling represents one of the most interesting universal 
amplitude ratios, being related to the fourth derivative of the free energy. 

Recently, in |[|, a new interesting approach has been proposed to evaluate 
this quantity in the case of integrable QFT's. The idea is that for these 
theories one has direct access to the so called form factors from which the 
renormalized coupling can be computed. In |]J the method was tested in 
the case of the 2d Ising model. The authors found the remarkably precise 
estimate 

g\ = 14.6975(1) (1) 

(see below for the precise definition of gl). 

The aim of this paper is to test this result by using a completely different 
method. By combining transfer-matrix methods and the exact knowledge of 
several terms in the scaling function of the free energy of the model we are 
able to obtain a precision similar to that of 0. Our result is 

gl = 14.69735(3), (2) 

which is in substantial agreement with the estimate of |1[. Given the subtlety 
of the calculations involved in both approaches, our result represents a highly 
non-trivial test of both methods. In performing our analysis we employ the 
same techniques which were used in in the study of the 2d Ising model in 
a magnetic field. 

This paper is organized as follows: We begin in sect. 2 by collecting some 
definitions and elementary results which will be useful in the following. Sect. 3 
is devoted to a discussion of the transfer-matrix results (and of the techniques 
that we use to improve the performance of the method). In sect. 4 we obtain 
the first 4 terms of the scaling function for the fourth derivative of the free 
energy, which enters in the estimate of g^ and finally, in sect. 5, we discuss 
the fitting procedure that we used to extract the continuum-limit value from 
the data. To help the reader to follow our analysis, we have listed in tab.||] 
the output of our transfer-matrix analysis. 
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2 General setting 



We are interested in the 2d Ising model defined by the partition function 

(7i=±l 

where the field variable a n takes the values {±1}; n = (no, ni) labels the sites 
of a square lattice of size L Q and L\ in the two directions and (n, m) denotes 
nearest neighbour sites on the lattice. In our calculations with the transfer- 
matrix method we shall treat asymmetrically the two directions. We shall 
denote no as the "time" coordinate and ri\ as the "space" one. The number 
of sites of the lattice will be denoted by N = L Li. The critical value of j3 is 

= (3 C = - log (V2 + 1) = 0.4406868... 

In the following we shall be interested in the high-temperature phase of the 
model in which the Z2 symmetry is unbroken, i.e. in the region (3 < f3 c . It is 
useful to introduce the reduced temperature t defined as: 

As usual, we introduce the free-energy density F(t, h) and the magnetization 
per site M(t, h) defined as 

F(t, h) = 1 log(Z(f , h)l M(t, h) = (5) 

The standard definition of the four-point zero-momentum renormalized cou- 
pling constant g 4 is 

i?(4) 

9^) = -372-' ( 6 ) 

A <52nd 

where x an d -^ 1( ' 4 ' ) are the second- and fourth-order derivatives of the free- 
energy density F(h,t) at h = 0: 



x(t)= W) 



F «> (t) = 



h=0 



{dhy 



(7) 

/i=0 



2 



and ^2nd denote the second moment correlation length. The second moment 
correlation length is defined by 

where d is the dimension (here d = 2) and 

fii = lim lim ^-V(m-n)'< a m a n > c . (9) 

The connected part of the correlation function is given by 

< c m cr n > c =< <y m o n > — < a m X a n > . (10) 
In particular, we are interested in the continuum-limit value g\ defined 

gl = hmg,(t). (11) 



as 



For t — > we have 



(12) 



from which it follows 



^4 = • ( 13 ) 



The amplitude A x is known exactly (see e.g. ref. 0): A x = 0.9625817322... 
The amplitude A^ 2 nd can also be computed exactly. Indeed, consider the 
exponential correlation length £ (inverse mass gap). For £ — > 0, it behaves as 
[|] A^~\ with A e = l/(4/3 c ) = 0.56729632855.... Using then § A 5 /A fi2n d = 
1.000402074..., we obtain finally A 5>2nd = 0.5670683251.... 

Our goal in the remaining part of this paper is to give a numerical estimate 
of A F (a). 

3 Transfer-matrix results 

We may have direct numerical access to by looking at the h dependence 
of the magnetization at fixed t. Expand as follows: 

h = b 1 M + b 3 M 3 --- (14) 
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we immediately see that 

h = Vx, &s = -^r, (15) 

so that 

F (4) = -663/61- (16) 

3.1 The transfer-matrix technique 

As a first step we computed the magnetization M of a system with L = 00 
and finite L\. The magnetization of this system is given by 

M = v%Mv , (17) 

where v is the eigenvector of the transfer matrix with the largest eigenvalue 
and M is a diagonal matrix with being equal to the magnetization of the 
time-slice configuration i. For a detailed discussion of the transfer-matrix 
method see e.g. ref. ||. 

We computed vq using the most trivial iterative method, 

= M\' (18) 

starting from a vector with all entries being equal. 

An important ingredient in the calculation is the fact that the transfer 
matrix can be written as the product of sparse matrices (see e.g. ref. |§). 
This allows us to reach L\ = 24 on a workstation. The major limitation is 
the memory requirement. We have to store two vectors of size 2 Ll . Since we 
performed our calculation in double precision, this means that 268 MB are 
needed. Slightly larger L\ could be reached by using a super-computer with 
larger memory space. 

For the parameters (3 and h that we studied, n < 200 was sufficient to 
converge within double-precision accuracy. 

3.2 The equation of state 

In order to obtain high-precision estimates of it turns out to be important 
to consider the external field h as a function of the magnetization rather than 
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the opposite. The advantage of the series ( |T4"D is that the coefficients — at 
least those we can compute — are all positive, and therefore, truncation 
errors are less severe than in the case of m(h). 

There is no sharp optimum in the truncation order. After a few numerical 
experiments we decided to keep in eq. fli~4]) the terms up to 615 M 15 : 

h(M) = b x M + 63 M 3 + • • • + &15M 15 . (19) 

In order to compute the coefficients 61, 63, 615 we solved the system of 
linear equations that results from inserting 8 numerically calculated values 
of the magnetization M(h\), M(/i2), M(h%) into the truncated equation 
of state ([19]). Here we have chosen hj = j hi. 

The errors introduced by the truncation of the series decrease as hi de- 
creases, while the errors from numerical rounding increase as hi decreases. 
Therefore, we varied hi to find the optimal choice. For a given value of (3 we 
performed this search only for one lattice size L\. (Typically Li = 18). From 
the variation of the result with hi we can infer the precision of our estimates 
of hi. For example, for (3 = 0.37, we get bi with 14 significant digits and 63 
with 12 significant digits. 

3.3 Extrapolation to the thermodynamic limit 

From the transfer matrix formalism it follows that for periodic boundary 
conditions and j3 7^ /3 C , the free energy density approaches its thermodynamic 
limit value exponentially in L\. Hence, also derivatives of the free energy 
density with respect to h and linear combinations of them should converge 
exponentially in Li to their thermodynamic limit value. Therefore, in the 
simplest case, one would extrapolate with an Ansatz 

b(Li) = b(oo) +cexp(-xLx) , (20) 

where b(Li) is the quantity at the given lattice size Li and b(oo) the ther- 
modynamic limit of the quantity. In order to obtain numerical estimates for 
6(00), c and x we have inserted the numerical result of b for the three lattice 



sizes Li, L\ — 1 and Li — 2 into eq. 20. It turns out that using this simple 



extrapolation, still a dependence of the result for 6(00) on Li is visible. This 
indicates that, with our numerical precision, subleading exponential correc- 
tions have to be taken into account. For this purpose we have iterated the 
extrapolation discussed above. 
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Table 1: Extrapolation of 63 to the thermodynamic limit for (3 = 0.37. 
Iterative procedure. The numbers in the top row give the extrapolation 
level. For the discussion see the text. 








1 


2 


3 


4 


13 


0.0459057204193 










14 


0.0463456447150 










15 


0.0467262982921 


0.049170965 








16 


0.0470483839889 


0.048819648 








17 


0.0473162239288 


0.048638691 


0.048446477 






18 


0.0475358880370 


0.048537476 


0.048409004 






19 


0.0477140164478 


0.048477931 


0.048392846 


0.048380598 




20 


0.0478571162217 


0.048441711 


0.048385463 


0.048379249 




21 


0.0479711755247 


0.048419155 


0.048381922 


0.048378658 


0.048378196 


22 


0.0480614838402 


0.048404863 


0.048380148 


0.048378366 


0.048378082 


23 


0.0481325804063 


0.048395686 


0.048379222 


0.048378214 


0.048378046 


24 


0.0481882778282 


0.048389731 


0.048378721 


0.048378128 


0.048378019 
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The iteration starts with b^(Li) which are the quantities b that have 
been computed by the transfer matrix for the lattice size L\. A step of the 
iteration is given by solving the system of equations 

b^{L x - 2) = c exp(-x(L 1 - 2)) + b^ x \L x ) 
6 {i) (L! - 1) = c exp(-x(L 1 - 1)) + & (<+1) (Li) 

= c exp(-xLi) + . (21) 

with respect to b^ l+1 \Li), c and x. In table [I] we give as an example the 
extrapolation of 63 at /3 = 0.37. In the second column we give the results 
obtained for the given L\. The stability of the extrapolation with varying L\ 
increases up to the fourth iteration. Further iterations become numerically 
unstable. 

As final result we took 63 = 0.04837802(3) from the 4 th iteration. The 
error was estimated from the variation of the results with L\. As a con- 
sistency check, we also extracted the thermodynamic limit by fitting with 
multi-exponential Ansatze. We found consistent results. The relative accu- 
racy of bi in the thermodynamic limit was in general better than that of 
63. 

In the second column of table |2| we give our final results for — F^H 11 ^ 2 at 
the P values that we have studied. For a discussion of the following columns 
see section 5. 



4 Scaling function for 

In this Section we shall study the asymptotic behavior of F^(t) for t — > 
following Ref. |§. With respect to |§, we have added the contributions due 
to the irrelevant operators. Here, we shall use the knowledge of the operator 
content of the theory at the critical point which can be obtained by using 
the methods of 2d conformal field theories. 

General renormalization-group (RG) arguments indicate that the free en- 
ergy of the model can be written as 

F(t,h) = f l(M0 + K |^ m! (^,{^}) 

+ l^logN/^^,{^}). (22) 
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Table 2: Extrapolation of ap 4 



p 






. 7 ,o 75 


ext. t 6 -' b 


ext. t 4 


0.200 


3.21111498292(2) 


4.202506 


4.358692 






0.250 


3.4721394466(1) 


4.286829 


4.369154 


4.383609 


4.382846 


0.280 


3.6225400346(6) 


4.321965 


4.373381 


4.381365 


4.380954 


0.300 


3.720514859(3) 


4.339751 


4.375424 


4.380395 


4.380141 


0.310 


3.768883189(7) 


4.347109 


4.376235 


4.380086 


4.379911 


0.320 


3.81687386(2) 


4.353517 


4.376917 


4.379774 


4.379635 


0.330 


3.86451569(5) 


4.359033 


4.377479 


4.379545 


4.379438 


0.340 


3.91183946(6) 


4.363716 


4.377934 


4.379382 


4.379302 


0.350 


3.9588780(3) 


4.367629 


4.378292 


4.379270 


4.379213 


0.355 


3.9823015(7) 


4.369315 


4.378439 


4.379239 


4.379196 


0.360 


4.0056667(8) 


4.370833 


4.378566 


4.379203 


4.379168 


0.365 


4.028978(1) 


4.372189 


4.378674 


4.379174 


4.379146 


0.370 


4.052240(2) 


4.373391 


4.378765 


4.379148 


4.379126 



Here Ff,(t, h) is a regular function of t and h 2 , the so-called bulk contribu- 
tion, ut, Uh, {uj} are the non-linear scaling fields associated respectively to 
the temperature, the magnetic field and the irrelevant operators, and yt, 
yh, {yj} are the corresponding dimensions. For the Ising model y t = 1, 
yh = 15/8. Notice the presence of the logarithmic term, that is related to 
a "resonance" between the thermal and the identity operator^. The scaling 
fields are analytic functions of t and h that respect the Z 2 parity of t and h. 
Let us write the Taylor expansion for Uh and u t , keeping only those terms 
that are needed for our analysis (we use the notations of 0): 

u h = h [1 + c h t + d h t 2 + e h h 2 + f h t 3 + 0(t\th% (23) 

u t = t + b t h 2 + c t t 2 + d t t 3 + e t th 2 + g t h 4 + f t t 4 +0{t 5 ,t 2 h 2 ). (24) 
Let us first discuss the contributions of the irrelevant operators. In generic 

In principle, logarithmic terms may also arise from additional resonances due to the 
fact that yj are integers or differ by integers from y^. They will not be considered here 
since these contributions either are subleading with respect to those we are interested in 
or have a form that is already included. 
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models their dimensions are usually unknown. In the present case instead, we 
may identify the irrelevant operators with the secondary fields obtained from 
the exact solution of the model at the critical point and use the corresponding 
RG exponents as input of our analysis. We shall discuss this issue in full detail 
in a forthcoming publication, let us only summarize here the main results of 
this analysis. It turns out that, discarding corrections of order 0(t 5 ), we 
have only two possible contributions: 

• The first one is due to terms TT , T 2 and T 2 (where T denotes the 
energy-momentum tensor). These terms would give a correction pro- 
portional to t 2 in the scaling function. 

• The second contribution is due to the L^ 3 L_ 3 I field from the Identity 
family and to L_ i e ) L_ 4 e from the energy family (where the L_^s are 
the generators of the Virasoro algebra). These terms give a correction 
proportional to t 4 in the scaling function. 

However, it turns out (see for instance the remarks of [0, |IU|) that in the 
infinite-volume free energy of the 2d Ising model the TT , T 2 and T 2 terms 
are actually absent^. Thus, from the above analysis we see that the first 
correction due to the irrelevant fields appears only at order t A . Therefore, 
since Uj/\ut\ v ^ Vt vanishes for t — ► 0, we can expand 

/ sing (x, { Zj }) = Y + (x) + u (t, h)u$X + (x) + OK 5 ), (25) 

where u (t, h) is an analytic function of t and h, and Y + , X + are appropriate 
scaling functions. The same expansion holds for / S i ng with different functions 
Y + , X + . Additional constraints can be obtained using the exactly known 
results for the free energy, the magnetization and the susceptibility in zero 
field. Since all numerical data indicate that all zero-momentum correlation 
functions diverge as a power of t without logarithms for t — > 0, we assume as 
in Ref. that Y + (x) is constant, i.e. Y+(x) = Y . The exact results for the 
free energy and the magnetization give then 

Pc , 23/? c 2 191/5 C 3 
Ch ~ 71' "16"' A "W2' (26) 

2 This conjecture is verified by the free energy and by the susceptibility at ft, = 
]10| and by the free energy F(0,h) Note that this is expected to be true only in 
the thermodynamic limit. In the finite-size scaling limit corrections that vanish like L^ 2 
are indeed observed pJ| . It is also not true for other observables, for instance, for the 
correlation length £. 
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c- Pc d- 7 ^ f- (27) 

ft ~^ (7) 

where we have adapted the numbers of to our normalizations, and Y = 
— 4/5^/7r. By making use of the expansion of the susceptibility, we obtain 
further 

lf>(0) = ^, * = -s|. < 28 > 

where Dq is the coefficient of the contribution proportional to t log \t\ in the 
susceptibility Numerically D = 0.04032550 . . ., so that b t = -0.0407708 . . . 
Nickel [10[] has also conjectured, on the basis of the numerical analysis of the 
high-temperature series of the susceptibility, that e t = b t j3 c \f2. 

Using the results presented above, and taking four derivatives of the free 
energy we obtain 

f& = r n / 2 (a F4 (i) +tV 4 (t)iog|t|) + r u ' 4 (b F i(t) +t A b F4 (t)log\t\) 

+c F4 (t) +c F4 (*)log|t|, (29) 

where ap^t), b F <±(t), c F ^{t), a F ±(t), b F ±{t), and c F ±(t) are analytic functions. 
Using Eqs. fl26|) and (^7]), we can compute the first terms in the Taylor 
expansion of a Fi (t). By direct evaluation we find 

= y^-^Me^yo^. ,30, 
From Eq. (|30|) , we immediately identify 

r| 4) (o) = A F(4) . (3i) 

Analogously, a direct calculation shows that 

b F4 (0) = -21 bt F| 2) (0) = 0.8241504.... (32) 

The contributions proportional to c F 4(t) and c F ±(t) give corrections of order 
t 11 / 2 which will be neglected in the following. 
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Putting together the various terms, we end up with the following expres- 
sion for the scaling function: 

F^t 11 ' 2 = A FW (l+p 1 t + p 2 t 2 + p 3 t 3 ) 

+p 4 t 11/4 + p^ 4 + p 6 t 4 + p e t 4 log \t\ + p 7 t 19 ^ 4 + 0(t 5 ) (33) 

where 

33 

pi = ^ = -0.46741893... (34) 

p 2 = — — = 0.052597147... (35) 
48 

29 3 3 

Ps = ^= = 0.054843243... (36) 
p 4 = -21b t Yi 2 \0) = 0.8241504.... (37) 

and P5,Pe,p6,P7 and A F ^ are undetermined constants which we shall try to 
fix in the next section. 



5 Analysis of the data 

The aim of this section is to obtain a numerical estimate for A F (a) by fitting 
the data reported in tab .0 with the scaling function (|3"3"|). The major problem 
in doing this is to estimate the systematic errors involved in the truncation 
of the scaling function. To this end we performed two different types of 
analysis. Let us see in detail the procedure that we followed. 



5.1 First level of analysis 

We first performed a rather naive analysis of the data. In table ^| we include 
step by step the information that we have gained in the previous section. In 

n/2 

the third column of table we have multiplied — F^ A > by , fe ^ 4 , where Uh 



and ut are given by eqs. (|23|j24]) . We see that the variation from (3 = 0.30 
to 3 = 0.37 of the numbers in column three is reduced by a factor of about 
10 compared with column two. In column four we add &F4(0) t 11 ^ 4 to the 
numbers of column three. Again we see that the variation of the numbers 
with 3 is drastically reduced in column four compared with column three. 
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Since we do not know the coefficients of higher order corrections exactly 
we have to extract them from the data. In the last two columns of table |2] 
we have extrapolated linearly in t x , with x = 3.75 in column 5 and x = 4 
in column 6. For the extrapolation we used neighboring /3-values (e.g. the 
value quoted for (3 = 0.37 is obtained from the extrapolation of the data for 
P = 0.365 and (3 = 0.37). 

We see that the result of the extrapolation does not vary very much when 
the exponent is changed from 15/4 to 4. Also the numbers given in column 
5 and 6 are much more stable than those of column 4. 

From this naive analysis we conclude that of 4 (0) = 4.3791(1), where the 
error bar is roughly estimated from an extrapolation of column 5 with £ 4 . 

In the next section we shall try to include the higher order corrections in 
a more sophisticated fitting procedure. 

5.2 Second level of analysis 

We made three types of fits: 

fl] In the first we kept Apw, p 5 and p 6 as free parameters. 

f2] In the second we kept A F n), p§, p$ and p-j as free parameters. 

f3] In the third we kept A F (a), p 5 , p 6 and p 6 as free parameters. 

These are the only choices allowed by the data. If we neglect also p§ we 
can never obtain an acceptable confidence level (in fact we know that pe is 
certainly different from zero and our data are too precise to allow such an 
approximation). If we add further terms, like a power of t 5 for instance, or 
try to fit simultaneously P5, pe, p% and p-j it always happens that some of the 
amplitudes are smaller than their statistical uncertainty signalling that our 
data are not precise enough to allow for five free parameters. 

In order to estimate the systematic errors involved in the estimate of A F (±) 
we performed for all the fitting functions several independent fits trying first 
to fit all the existing data (those listed in tab.|2|) and then eliminating the 
data one by one, starting from the farthest from the critical point. Among 
the set of estimates of the critical amplitudes we selected only those fulfilling 
the following requirements: 
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1] The reduced x 2 of the fit must be of order unity. In order to fix precisely 
a threshold we required the fit to have a confidence level larger than 
30%. 

2] For all the subleading terms included in the fitting function, the am- 
plitude estimated from the fit must be larger than the corresponding 
error, otherwise the term is eliminated from the fit. It is exactly this 
constraint which forbids us to take into account fits with more than 
four free parameters. 

3] The amplitude of the n th subleading field must be such that when it is 
multiplied for the corresponding power of t, (for the largest value of t 
involved in the fit) it gives a contribution smaller than that of the (n — 
l) th subleading term. This is intended to avoid artificial cancellations 
between subleading fields. 

Among all the estimates of the critical amplitude A F w fulfilling these re- 
quirements we select the smallest and the largest ones as lower and upper 
bounds. 

The results of the fits are reported in tab.Q, f| and |. We report all the 
combinations of input data which fulfill requirements [l]-[3] . In the tables 
we also report the best fit value of p$. All the fits were performed using the 
double-precision NAG routine G02DAF. 



A F (4) 


P5 


d.o.f. 


-4.3791092(9) 


0.236(3) 


4 


-4.3791065(9) 


0.220(12) 


3 


-4.379095(2) 


0.142(18) 


2 



Table 3: Fits of type [fl] fulfilling the requirements 1-3. In the first column 
the best fit results for the critical amplitude (with the error induced by the 
systematic errors of the input data in parenthesis), in the second column the 
best fit value ofps, in the last column the number of degrees of freedom (i.e. 
the number of data used in the fit minus the number of free parameters) . 
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A F (4) 


P5 


d.o.f. 


-4.3791003(1) 


0.0486(7) 


7 


-4.3791001(3) 


0.047(3) 


6 


-4.3791006(8) 


0.053(9) 


5 


-4.3791022(14) 


0.079(22) 


4 



Table 4: Same as tab.f|, but with fits of type [f2]. 



A F (4) 


P5 


d.o.f. 


-4.3790944(1) 


-0.692(1) 


8 


-4.3790942(3) 


-0.696(7) 


7 


-4.3790961(6) 


-0.634(17) 


6 


-4.3790980(10) 


-0.56(4) 


5 


-4.3791005(16) 


-0.41(8) 


4 



Table 5: Same as tab.^, but with fits of type [f3]. 



Looking at the three tables and selecting the lowest and highest values of 
A F (4) we obtain the bounds 



- 4.379093 > A FW > -4.379110 



from which, using eq . (|T3|) , we obtain 

gl = 14.69735(3) 



(38) 



(39) 



which we consider as our best estimate for g\. As anticipated in the in- 
troduction, this result is in substantial agreement with the estimate of [p]]. 
Notice however that the error quoted in eq.(|3"9"|) should not be considered as 
a standard deviation. It rather encodes in a compact notation the systematic 
uncertainty of our fitting procedure. 

We can compare the estimate ([39]) with previous numerical determina- 
tions. The analysis of high-temperature expansions gives g\ = 14.694(2), 
Ref. PJ and gl = 14.693(4), Ref. p[ while Monte Carlo simulations give 
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g% = 14.3(1.0), Ref. [[TJ, and g\ = 14.69(2), Ref. §. These results agree 
with our estimate (|39|) , which is however much more precise. 

It is clear from the data (see the second column of tab. [3], f| and ^) that 
the uncertainty on A F (a) is mostly due to the fluctuation of p 5 . If one would 
be able to fix exactly also p§, the precision in the determination of g± could 
be significantly enhanced. 
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